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5_i ' Abstract 



We examine the potential of the COMPASS experiment at CERN to study color 
transparency via exclusive vector meson production in hard muon-nucleus scattering. It 
is demonstrated that COMPASS has high sensitivity to test this important prediction 
of perturbative QCD. 



1 Expanded version of the talk presented at the Workshop on "Nucleon Structure and Meson Spectroscopy" 
physics to be studied at the COMPASS experiment at CERN, Dubna, Russia, 10-11 October 2000. The 
conference www site is http://wwwcompass.cern.ch/compass/workshops/dubna2000/index.html. 
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1 Introduction 



Color transparency (CT) is a phenomenon of perturbative QCD (pQCD), whose char- 
acteristic feature is that small color-singlet objects interact with hadrons with small cross 
sections [1, 2, 3, 4]. Cross section for the interaction of such small object, or small size 
configuration (SSC), with a hadron target has been calculated in QCD using the factor- 
ization theorem [5, 4, 6, 7]. These QCD calculations confirmed the hypothesis of F. Low 
[1] on smallness of the cross section for the interaction of SSC with a hadron, if the gluon 
density in the hadron is not very high (moderately small x). They also predict a related 
phenomenon of color opacity when the gluon density becomes very large (at small x) and 
SSC interacts with the hadron with large cross section [2, 4]. 

The prerequisite for observing CT is to select a sample containing SSC's via a hard 
process (i.e. with large Q 2 , high p t , or large produced mass). To suppress non-perturbative 
(not SSC) background different additional restrictions should be imposed depending on 
the process. For instance for hard exclusive p° leptoproduction, in addition to large Q 2 , 
selection of the longitudinally polarized mesons is required. To 'measure' the SSC-nucleon 
cross sections one should study absorption of the SSC propagating through nuclear matter. 
In order to clearly observe CT it is necessary that the SSC lives long enough to traverse 
distances larger than the size of the target nucleus. Another requirement is that SSC stays 
small while propagating through the nucleus. These requirements, which are quantified in 
terms of the coherence length and the formation length, are discussed later. 

Various processes have been proposed to study CT phenomenon. 

(a) Coherent vector meson (J/0, p, (ft) production at small t. For near-forward (t « 0) 
coherent production and complete CT, one expects that the cross section for the 
nuclear target is related to that for the nucleon target by 

^ (7 M -> VA) U = A 2 ^(YN VN) U . (1) 

Gluon shadowing/antishadowing in nuclei is neglected in this formula, but for mod- 
erately small x(> 0.01) this effect is expected to be small. Usually in experiments the 
i-integrated coherent cross sections are measured, for which CT predicts an approx- 
imate A 4 / 3 dependence. Using more realistic nuclear form factors one predicts A 1A0 
[8]. 

Otherwise, if the color-singlet objects, while propagating in nuclear matter, interact 
with nucleons with large cross sections, the expected yl-dependence is weaker; A 2 / 3 
for the cross sections comparable or larger than the pion-nucleon cross section. 

(b) Incoherent (quasi-elastic) vector meson (J/ift, p, (ft) production on nuclei. For complete 
CT and neglecting effects of the gluon shadowing/antishadowing in the nuclei one 
expects 

^( 7 M -> VN{A - 1)) = A^tfN VN) . (2) 

Here, for the nuclear target the meson V is produced on a single nucleon of the nucleus 
and (A — 1) denotes the system of spectator nucleons. 

(c) Coherent or incoherent production of excited vector meson states ift' or p'. CT predicts 
the same A dependence of J /ift and ift' , or p and p' . This contrasts with the naive 
expectations whereby one may expect larger absorption for excited mesons since they 
are larger. 
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(d) Coherent diffractive dissociation of hadrons or photons into high p t di-jets. Such 
process probes the small transverse-size component of the projectile wave- function 
as well as CT effects. For CT the t-integrated cross section of coherent diffractive 
production of high pt di-jets has the same ^4-dependence as for the processes (a). 
Using more realistic wave functions one predicts for asymptotically high energies A a 
dependence, where a is in the range 1.45 - 1.60, depending on p t [4]. 

(e) Coherent vector meson production on light nuclei (deuteron, helium) in the large 
t range, where the effects of double (multiple) scattering are important. CT will 
suppress the double (multiple) scattering contribution to the differential cross section 
da/dt [4, 9]. 

(f) ^4-dependence of the fraction of large rapidity gap (diffractive) events. CT predicts 
an A independent fraction; otherwise it will grow with A [10]. 

(g) Large Q 2 quasielastic (e,e'p) scattering on nuclei, eA — > ep(A — 1). In the CT limit 
the cross section will be proportional to A, as for process (b) [11]. 

(h) Large t quasielastic (p, 2p) scattering on nuclei, pA — > pp(A-l). Also for this process, 
for complete CT the expected A dependence is like for process (b) [11]. 

In searches for CT a commonly used quantity, measured in experiments, is the nuclear 
transparency 

T = £-, (3) 

A(J N 

which is the ratio of the cross section per nucleon for a selected process on a nucleus A 
to the corresponding cross section on a free nucleon. For the incoherent processes CT 
predicts T ~ 1 independently of A. For the large absorption in the nuclear matter T will 
be smaller than unity and will decrease with A. Although the nuclear transparency could 
be defined also for coherent processes, usually cross sections for coherent production on 
different nuclei are compared directly. As mentioned above, for i-integrated coherent cross 
sections CT predicts a A oc A 4 / 3 , whereas for a larger absorption the ^-dependence is weaker. 

In the following discussion we concentrate on exclusive vector meson production 
(VMP), the processes (a) - (c), which could be studied at the COMPASS experiment [12, 13]. 
Exclusive meson production on the nucleon, free or bound in the nucleus, can be viewed 
as proceeding according to the diagram shown in Fig. 1. Kinematic variables used in this 
paper are listed in Table 1. 

One of predictions of pQCD is that at high Q 2 the longitudinally polarized virtual 
photons 7£ fluctuate into hadronic components, e.g. qq pairs, whose transverse size b = 
| f± q — fj_q | decreases with Q 2 , b oc (Q 2 ) -1 / 2 . At large Q 2 the values of b are significantly 
smaller than the size of the nucleon. For instance, at Q 2 = 10 GeV 2 and at x = 0.01 
b U u — b d g ~ 0.3 fm [6], to be compared with the size of p meson b p ~ 1.4 fm. For 
transversely polarized photons, in addition to the small size fluctuations, non-perturbative 
large size components may be expected, even at reasonable large Q 2 (~ 10 GeV 2 ). 

Another pQCD prediction is that the heavy quarks fluctuations, e.g. cc, of the virtual 
photon have small transverse size already for quasi-real production; b cS = 0.22 fm at Q 2 ~ 0. 

The total cross section for the interaction of a small size qq pair with the nucleon is 
given in pQCD by the formula [5] 



7T 2 

<7qq, N = ^-b 2 a s {Q 2 ) x g(x, Q 2 ) , (4) 
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Table 1: Kinematic variables used in the text. 



k 


fAlir TT"1 OlTVl OT1 T 11 TV1 t~~iT TrlO 1 VI /**1 f\ OV1T miir~1Y1 

1U LL1 IllUIllCllL Lllll UL LUG IIICILLCIIL I11LIU11, 


K 


four-momentum of the scattered muon, 


P 


four-momentum of the target nucleon, 


V 


fmir-mOTMpntiiTn nf trip vprtnr Tnpsnn V 


q — k - k' 


four-momentum of the virtual photon, 


-Q 2 = q 2 


invariant mass squared of the virtual photon *y* , 


v = (n ■ a) /Mr, 


energy of the virtual photon in the laboratory system, 
Mp is the proton mass, 


x = Q 2 /{2M p v) 


Bjorken scaling variable, 


v = {p- q)/(p ■ k) 

W 2 = (p + q) 2 
t = (q-v) 2 


fraction of the lepton energy lost in the laboratory system, 


total energy squared in the 7*— TV system, 


four-momentum transfer to the target, 


p\ 


transverse-momentum squared of the vector meson 
with respect to the virtual photon direction, 


m v = (f 2 )' 


invariant mass of the vector meson V, 


Ml = (p + q-v) 2 
I = {Ml-M^)IW 2 


missing-mass squared of the undetected recoiled system, 


inelasticity, 


z = (p- v)/{p ■ q) 


fraction of the virtual photon energy in the laboratory sy- 
stem taken by meson V , 



where g(x, Q 2 ) is the gluon distribution function in the nucleon. The cross section for the 
interaction of qq pair of b = 0.3 fm with the nucleon is about 3 mb at x = 10~ 2 [6]. At 
very small x the gluon distribution g in the nucleon increases, which leads to an increase 
of the cross section a q q,N- For instance, for b = 0.3 fm cr q q,N is about 18 mb at x = 10~ 5 . 
Therefore, even for the small qq objects the cross section becomes large at sufficiently small 
x. This phenomenon is mentioned earlier color opacity and it could be studied at small x 
at future electron-nucleus colliders. 

In order to study CT for the exclusive production of light quark mesons (like p or p') 
one should select large Q 2 , moderately small x and longitudinally polarized mesons. For 
similar studies for the production of J/tp or ip 1 one should just select moderately small x. 

In studies of CT an important role is played by the coherence length and the formation 
length. The coherence length l c (sometimes referred to as propagation or interaction length) 
is defined as the distance traversed by the qq fluctuation of the virtual photon in the target 
nucleon/nucleus system and is given by 

c Q 2 + M 2 M N x' 1 ' 

Here (3 = Q 2 /{Q 2 + M 2 ), M is the invariant mass of the qq fluctuation, and M N is the 
mass of the nucleon. For l c values smaller than the size of the nucleus the life-time of the qq 
fluctuation is short and a chance for the hadronic fluctuation to interact in nuclear matter 
decreases. In particular this significantly inhibits coherent production. For incoherent 
production the effects due to the small l c values mimic the Q 2 dependence of the nuclear 
transparency predicted by CT [14]. 

The formation length is the distance in the target system needed for the qq fluctuation, 
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which scattered on a nucleon, to develop into a hadron h. It is equal to 

h = — V -s— > ( 6 ) 

J m ft Am 

where m h is the mass of the hadron h and Am is the mass difference between the hadron h 
and its lowest orbital excitation. For small If values the qq fluctuation evolves quickly into 
a full-size final hadron and the absorption of the final hadron h in the nucleus plays a role. 

Experimentally the effects of the coherence length and the formation length were 
analysed in detail by the HERMES experiment for incoherent exclusive production of p° on 
nuclei [15] and were shown to play an important role at the small l c values. 

Therefore, for a clean demonstration of CT effects the optimal conditions are when 
the values of l c and If exceed the size of the target nucleus. If the above is not feasible, the 
variation of nuclear transparency T(A, Q 2 ) with Q 2 should be studied for different, fixed 
values of coherence length. This way, a change of T(A) between low and high Q 2 values 
could be associated with the onset of CT, and not with varying l c . 

The t-dependence of the cross section for exclusive meson production on the nucleon is 
approximately exponential, e h ^~ tmin \ where t m i n ~ —M 2 x 2 (l+M 2 /Q 2 ) 2 and M is the mass 
of the produced meson. For soft QCD processes the slope parameter b should be independent 
of Q 2 , depending only on the energy. This results from the Regge model [16]. For hard 
VMP, if SSC's are important, one predicts [17] that the slope b will significantly decrease 
with increasing Q 2 approaching a universal value related to the nucleon radius. This is so 
because the transverse interquark distances in the SSC decrease with increasing Q 2 . This 
prediction agrees with high energy data on hard diffractive vector meson production on the 
proton [18, 19, 20]. 

For VMP on a nucleus A the i-dependence of the cross section is approximately 
reproduced by a sum of two exponential functions. The peak at the lowest t values, with 
the slope proportional to the nucleus squared radius <R\>, is mostly due to the coherent 
production, whereas at somewhat larger t the incoherent production on quasi-free nucleons 
dominates and the slope is equal to that for the production on free nucleons. Applying cuts 
on t allows to select samples of events, which are strongly dominated by either coherent or 
incoherent production. It was successfully demonstrated in Refs [8, 23, 26]. 

In order to observe maximal coherent exclusive meson production for which the whole 
nucleus contributes, it is necessary to satisfy the condition 

Kin <i?l>/3| « 1. (7) 

In the COMPASS experiment for low x values in the range 0.006 - 0.02 and for Q 2 ~ 
2 — 10 GeV 2 , one has \t m i n < R\ > /3 | ~ 10x 2 A 2 / 3 and the condition (7) is satisfied even 
for lead (A = 207). 

Experimental searches for CT started more than a decade ago and encompass various 
processes: large t quasielastic (p, 2p) scattering [21], large Q 2 quasielastic (e, e'p) scattering 
[22], J/0 photoproduction [8] and J/tp muoproduction [23], exclusive p° leptoproduction 
[24, 25, 15] and coherent diffractive dissociation of the pion into two high-p t jets [26]. 

The pioneering studies of CT in large t quasi-elastic proton scattering [21] found a rise 
of the nuclear transparency as the beam energy increased from 6 to 9 GeV and a decrease 
at higher energies. The large Q 2 quasi-elastic electron scattering studies [22] did not show 
Q 2 dependence. The explanation of these results in terms of CT is still debatable. 

Strong recent evidence for CT comes from Fermilab E791 experiment on the A-de- 
pendence of coherent diffractive dissociation of pions into two high-p t jets [26]. Also the 
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E691 results on ^4-dependence of coherent J/tp photoproduction [8] and the NMC measure- 
ments of ^4-dependence of coherent and quasielastic J/if) muoproduction are consistent with 
CT. Measurements of the nuclear transparency for incoherent exclusive p° production by 
Fermilab experiment E665 give a hint for CT. However, due to the low statistics of that 
data at high Q 2 , it was not possible to disentangle effects of decreasing l c at high Q 2 and 
to demonstrate CT unambiguously. 

2 Experimental method 

We propose to study CT via exclusive vector meson production pA — > p VA (coherent) 
and pA — > p VN(A — 1) (incoherent) on various nuclei A and optionally also on a proton 
or deuteron target. As a primary objective we propose to study the production of the 
following mesons V: p°, J/ip, <p, ip\ip{2S)) and / o'(p(1450),p(1700)). Also investigations of 
the production of other mesons will be possible. The preferable decay modes are those into 
the charged particles: p° — > ir + ir~, J/0 — > e + e~ / p + p~, 4> — > K + K~, tp' — > J/0 7r + 7r~, p' — > 
7T + 7r + 7r _ 7r _ . Two or more nuclear targets will be used. An additional proton (deuteron) 
target would be beneficial. 

Our proposal is to complement the initial setup of the COMPASS [12, 13] for the muon 
run with the polarized target, by adding two thin nuclear targets of lead and carbon of 17.6 
g/cm 2 each. The carbon target will be a cylinder 8 cm long and of 3 cm of diameter. The 
lead target will consist of 4 discs of 3 cm of diameter, distributed over length of 8 cm. Only 
one nuclear target will be exposed to the muon beam at a time, with frequent exchanges 
(every few hours) of different targets. The nuclear targets will be located downstream of 
the polarized target, at the end of the solenoid magnet and before the first tracking detector 
(first micromega chamber). 

The high-intensity high-energy incident muon-beam will impinge on the polarized 
target and a downstream thin nuclear target. The momenta of the scattered muon and of 
the produced charged particles will be reconstructed in the two magnetic spectrometers, 
using the magnets SMI and SM2, instrumented with micromega chambers, drift chambers, 
GEM detectors, straw chambers, multiwire proportional chambers and scintillating fibers. 
Adding a recoil detector which will surround the nuclear target and register slow particles 
emitted from it may not be possible for the present setup due to the limited space. However, 
it would be advantageous and possible for dedicated runs taken with a modified setup. In 
the following we assume that the recoil particle(s) remain(s) undetected, and in order to 
select exclusive events one has to rely only on the kinematics of the scattered muon and the 
produced meson. 

The trigger will use the information from hodoscopes registering the scattered muons, 
and from calorimeters registering deposits of energy of the particles in the final state. For 
VMP reactions triggering only on the scattered muon is in principle possible. However, for 
better efficiency of the trigger, especially at small Q 2 , it will be useful for certain processes 
to require in addition a minimum energy deposit in the calorimeters. 

The off-line selection of exclusive events for the production of different mesons will 
be similar to that described for p° production in Ref. [25, 27, 28]. In particular the 
discrimination of non-exclusive events will be done by applying cuts on the inelasticity / 
(for the definition see Table 1). In Fig. 2 the inelasticity distribution is shown for the SMC 
p° sample [29] for the events with the invariant mass in the central part of the p° peak. 
For the inelasticity distribution the peak at / = is the signal of exclusive p° production. 
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Non-exclusive events, where in addition to detected fast hadrons, slow undetected hadrons 
were produced, appear at / > 0. For the cut —0.05 < I < 0.05 defining the exclusive 
sample the amount of the residual non-exclusive background for the SMC experiment was 
up to about 10% at large Q 2 . The kinematical smearing in / and the width of the elastic 
peak in COMPASS is expected to be about the same (cf. Sect. 3.2) as that shown in 
Fig. 2 for the SMC experiment. Although the smearing will be similar, we expect the level 
of the non-exclusive background to be lower in COMPASS due to the wider angular and 
momentum acceptance coverage for final state hadrons. In addition, with larger statistic 
in COMPASS it will be possible to apply more tight inelasticity cuts, further reducing the 
background. The effect of this residual background on various observables will be studied 
by varying the inelasticity cuts, similarly as was done in Ref. [30]. 

The selections of coherent or incoherent production will be done on a statistical basis, 
using the t-distribution; at the lowest 1 1 | values coherent events predominate, whereas at 
somewhat larger 1 1 | there is almost clean sample of incoherent events. 

Separation of the p° samples with the enhanced content of longitudinally or trans- 
versely polarized mesons will be done by applying cuts on the measured angular distributions 
of pions from the decays of the parent p°. 

The minimal covered Q 2 range is expected to be 0.05 < Q 2 < 10 GeV 2 . For the med- 
ium and large Q 2 values (Q 2 > 2 GeV 2 ) the range 0.006 < x < 0.1 will be covered with 
good acceptance. 

The basic observable for each process studied will be the ratio of the nuclear trans- 
parencies for lead and carbon, R T = T^/Tq = (crpb/^4pb)/(ccMc)- Due to the proposed 
geometry of the targets, the acceptances will cancel in the ratio i? T . Also the absolute 
beam flux measurement will not be necessary for the ratio R T , provided that the relative 
determination of the beam fluxes for the exposures with different target materials could be 
done, e.g. by counting DIS events originating in the polarized target. The measured ratio 
R T should be corrected for different losses of events in lead and carbon, which are due to 
the secondary interactions in the targets. They will be estimated from the MC simulations. 

3 Simulation of exclusive p° events 

In this section we describe details of the simulation of exclusive coherent (pA — ► p p°A) 
and incoherent (pA — > p p°N(A — 1)) p° production in the COMPASS experiment with the 
carbon and lead targets. The simulations were done with a dedicated fast Monte Carlo 
program which generates deep inelastic exclusive p° events with subsequent decay p° — > 
7r + 7r~. At this stage there was no attempt to include any background in the event generators. 
Here we describe the event generator as well as the treatment of different experimental 
aspects: losses due to the secondary interactions of pions, propagation through the magnetic 
fields, angular and momentum resolutions, muon trigger acceptance, acceptance for final 
state pions and efficiency of tracks reconstruction. 

Radiative corrections for exclusive p° production are expected to be similar to those 
for the NMC experiment, which were in the range of 2% to 5% [31]. These corrections were 
not included in the present simulations. 
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3.1 Event generator 

First we present the used parameterization of the cross section for the production on 
the free nucleon, p,N — > pi p°N, with the subsequent decay p° — > 7r + 7r~: 

<7»N^ P °N = T T (Q 2 , V) • <T%^ N {Q 2 , ") • F (pI COS ^, 0) , (8) 

where 9 and </> are, respectively, the polar and azimuthal angles of ir + from the decay, 
calculated in the parent p° center-of-mass system, with respect to the direction of flight of 
p°, T T is the flux of transverse virtual photons 



a is the fine-structure constant, E^ the muon-beam energy and e is the virtual photon 
polarization given by 

1 - - 



Q 2 



e = 7 ^ n2 . (10) 



The virtual photon cross section was parametrized as 



1.96 



<,V vw (Q 2 ,,) = 27.4nb.(^j . (11) 

This is the NMC parametrisation of the cross sections per nucleon for exclusive p° produc- 
tion on carbon [25]. As the NMC data shows little A-dependence of the virtual photon (t-'m- 
tegrated) cross section per nucleon, we use the same parametrisation of ^*^_, p o^ f° r car ~ 
bon and proton targets. This parametrisation does not apply at small Q 2 values, namely 
at Q 2 < 1 GeV 2 

The function F comprises the p 2 distributions of produced p° and the angular distri- 
butions of pions coming from its decay 

F = a L ■ f L (p 2 ) ■ W L (cos 6, 4>) + a T - f T (p 2 ) ■ W T (cos 9, <P) . (12) 

Here a L = r^, a T = 1 — Tqq, Tqq = Tqq(Q 2 , v) is the p° density matrix element, which can be 
identified as a fraction of longitudinally polarized (helicity = 0) p° mesons, and the indices 
L and T refer to longitudinally and transversely polarized p 0, s, respectively. The fraction 
Tqq can be expressed [32] by the ratio R of the cross sections for exclusive production by 
longitudinal and transverse virtual photons. We use a parametrisation of R given by [33] 
which reproduce data on exclusive p° production in a wide range of Q 2 . 
The p 2 distributions are described by 

f l {p 2 ) = b l e~ h ^ , (13) 

where i = L or T, b L = 4.5 + 4 • (0.5/Q 2 ) GeV -2 and b T = 8.5 GeV" 2 These parameteri- 
zations weighted by the fractions of longitudinally and transversely polarized mesons allow 
to reproduce reasonably the values of the effective slope b for the exclusive p° production 
measured at HERA and in the fixed-target experiments in a wide range of Q 2 . 
The angular distributions of the pions from p 0, s decays are given by 

W$(cos0, # = i-|{(l-i?) + (3P t -l)cos 2 9} , (14) 
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where i = L or T, P L = and P T = 1. For the nuclear targets we assume the same W- 
distributions as for the proton as suggested by [25]. Note that the distributions fi(p 2 ), 
W i (cosO, <p) and F(p 2 , cosO, 4>) are normalized to unity. 

The invariant mass of two decay pions was generated using the relativistic p-wave 
Breit-Wigner shape for the p resonance [34]. 

We relate the differential cross sections for the proton 

= Tt ' ' ai ■ ' (15) 

to these for coherent and incoherent production on the nucleus A by 

(*£) r ^..<*>».(%.\, (16) 

respectively. Here < R\ > is the mean squared radius of the nucleus, A eS ^ coh and A cS ^ inc 
take account of nuclear screening for the coherent and incoherent processes, correspondingly. 
The cross section for incoherent exclusive meson production is summed over all final states of 
the recoiling system, i.e. it is given for the so called closure approximation. The suppression 
of the incoherent cross section at small t due to the Pauli blocking is neglected here. We 
used the approximation t — i min ~ ~Pt, where |t min | is the minimal kinematically allowed 
1 1 1 value for given W 2 , Q 2 , rriy and Mf . 

We generated the cross sections for two models. For the complete color transparency 
model (CT model) we used A eS coh = A eS inc = A. In another model we assumed a 
substantial nuclear absorption (NA model) and used A e ^ coh = A cS ^ inc = A - 75 . These could 
be compared to different predictions of the vector meson dominance model (VMD), which 
vary depending on the p meson-nucleon total cross section and on the coherence length l c 
(see e.g. [35]). The A dependences of the NA model are in the range of predictions of 
VMD, except for incoherent production at large l c , where VMD predicts stronger nuclear 
absorption. 



3.2 Simulation of the experimental effects 

The secondary hadronic interactions of the decay pions in the target were simulated. 
The assumed density of targets, p tgt , was 2.2 g/cm 3 . The interaction length Xf nt for pions 
in the target material, was assumed equal to 130 g/cm 2 for the carbon and 290 g/cm 2 for 
the lead target. For an exclusive p° event to be reconstructible, it was required that none 
of the decay pions underwent an inelastic hadronic interaction. 

The trajectories of charged particles were simulated taking into account the geometry 
of SMI and SM2 magnets of the COMPASS experiment. Homogeneous fields inside both 
magnets were assumed. For the SMI magnet the bending power jBdl was assumed equal 
to 1.0 T-m (independent of the beam energy), whereas for the SM2 magnet it was assumed 
equal to 2.3 T-m for a 100 GeV muon beam energy, and 5.2 T-m for a 190 GeV muon beam. 
After tracking of the produced charged particles through the detector, a flag was assigned 
to each particle telling how far it propagated in the COMPASS setup. 

Kinematic smearing of the beam, of the scattered muon and of the charged hadrons 
was simulated. Assumed values of dispersions of measured particle momenta and angles are 
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based on the experience of previous muon experiments at CERN, as well as on the results 
of studies at COMPASS. The relative error on the momentum, a(p)/p, was assumed equal 
to 0.5% for beam tracks, 0.75% for the tracks passing only the first magnet and 0.44% for 
the tracks passing the second magnet. The error on the angle of a particle was assumed to 
be 0.15 mrad. For a 190 GeV beam and the kinematic cuts listed in Section 3 the resulting 
smearing of the inelasticity / is about 0.018 and the smearing of the invariant mass of two 
pions is about 6 MeV. 

To simulate the trigger acceptance a trajectory of the scattered muon behind the 
second magnet was calculated, and the hits in the muon hodoscopes H4 and H5 were checked. 
Each of these two hodoscopes consists in fact of a few different hodoscopes (namely of the 
Ladders, the Primed System, the Unprimed System) but for the purpose of the present 
analysis we will consider just two cases, corresponding to different trigger acceptances. 
First, we assumed that only the Ladders and the Primed System are available. This trigger 
is called the Medium Q 2 range Trigger (MT). If in addition the Unprimed System is also 
implemented, the Full Q 2 range Trigger (FT) will result. The Q 2 dependence of the trigger 
efficiency, e tr , for these triggers was presented in [27]. For the MT trigger at 190 GeV beam 
e tr decreases quickly with Q 2 from 0.7 at Q 2 = 2 GeV 2 to about 0.1 at Q 2 = 10 GeV 2 The 
acceptance is several times smaller for this trigger at 100 GeV beam. For the FT trigger 
the Q 2 dependence is weaker and the trigger acceptance is higher; it is always bigger than 
0.5 for Q 2 < 70 GeV 2 at 190 GeV beam energy and for Q 2 < 20 GeV 2 at 100 GeV beam. In 
conclusion, for the MT trigger the data taking with the higher beam energy seems the only 
acceptable choice, whereas for the FT trigger running at both beam energies is feasible, 
although the covered Q 2 range is larger at the higher beam energy. 

For a p° meson to be accepted it was required that each pion from its decay was 
emitted in the laboratory at an angle within the acceptance of the SMI magnet, and that 
its momentum was bigger than 2 GeV. 

Based on the preliminary results of the track reconstruction by the programs devel- 
oped at COMPASS, simple and rather conservative assumptions were used to simulate the 
efficiency of tracks reconstruction. For the tracks seen only in the first spectrometer the 
single track reconstruction efficiency was assumed equal to 0.8 for the momentum range 
p > 2 GeV, and for the tracks observed also in the second spectrometer it was assumed 
equal to 0.95. The efficiency e rec to reconstruct all three tracks of the scattered muon and 
of two pions was assumed to be equal 0.7 • 0.95 3 + 0.3 • 0.95 2 • 0.80. This assumption was 
motivated by the observation that for 70% of accepted DIS exclusive p° events all three mea- 
sured tracks are seen in the second spectrometer, while for the remaining 30% events the 
scattered muon and the fast pion are seen in both first and second spectrometers, whereas 
the slow pion is observed in the first spectrometer only. 

4 Results on exclusive p° production 

Due to the higher trigger efficiency and larger Q 2 range at higher beam energy we 
considered 190 GeV muon beam. The simulations were done independently for the carbon 
(A = 12) and lead targets (A = 207), and for two triggers (MT and FT). For each target 
and each trigger we assumed two different models describing the nuclear effects for exclusive 
p° production: CT model and NA model (cf. Section 3.1). 
The kinematic range considered was the following: 

2 < Q 2 < 80 GeV 2 , (18) 
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Table 2: Numbers of accepted events for two considered models of the nuclear absorption. 



model 


N c 




CT 


70 000 


200 000 


NA 


28 000 


20 000 



35 < v < 170 GeV . (19) 

The upper cut on v was chosen to eliminate the kinematic region where the amount of 
radiative events is large, whereas the lower one to eliminate the region where the acceptance 
for pions from p° decay is low. 

The total efficiency e tot to observe exclusive p° events results from: the acceptance of 
the trigger (e tr ), acceptance to detect the pions (e had ), efficiency for tracks reconstruction 
(e rec ), cut on the invariant mass of two pions, 0.62 < M m < 0.92GeV 2 , used for the selection 
of the samples (e mass ), and efficiency for an event to survive the secondary interactions (e sec ). 
The contributions of all these effects to e tot are similar to those presented in [27] for the 
polarized target, except for the effects of the secondary interactions. The approximate value 
of e sec is equal to 0.87 for the carbon target and 0.94 for the lead target. The total efficiency 
e tot is about 0.48 for the FT trigger and about 0.30 for the MT trigger. 

The total expected cross section for exclusive p° production on the nucleon is 

v^uflN = L °»n^n(Q 2 i ") dQ 2 du , (20) 



f^N^pON ~ J u pN->w< 

min ^min 

where the kinematic range was defined before. The value of cr to ^^ is 283 pb. For 
nuclear targets the corresponding values depend on A and on the assumed model for nuclear 
absorption. 

The expected muon beam intensity will be about 10 8 /s during spills of length of about 
2 s, which will repeat every 14.4 s. With the proposed thin nuclear targets, each of 17.6 
g/cm 2 , the luminosity will be C = 12.6 pb -1 • day" 1 . The estimates of the numbers of 
accepted events were done for a period of data taking of 150 days (1 year), divided equally 
between two targets. An overall SPS and COMPASS efficiency of 25% was assumed. The 
numbers of accepted events for the carbon and lead targets, assuming the two models for 
the nuclear absorption mentioned earlier, are given in Table 2. 

The distributions of accepted exclusive events as a function of x, Q 2 and v are similar 
to those presented in [27]. 

In Fig. 3 we present the p 2 distributions for both targets. We observe clear coherent 
peaks at small p 2 (< 0.05 GeV 2 ) and less steep distributions for the incoherent events 
at larger p\. The arrows at the top histograms indicate the cut p\ > 0.1 GeV 2 , used 
to select the incoherent samples. The contribution from coherent events is negligible in 
these samples. For the middle and bottom histograms the arrows indicate the cut p 2 < 
0.02 GeV 2 , used to select the samples which are dominated by coherent events — the 
so called coherent samples. For the samples defined by the latter cut the fraction of the 
incoherent events is at the level of up to 10%, depending on the nucleus and on the model 
for nuclear absorption. We plan also to use the standard method to determine coherent and 
incoherent components, by fitting p 2 distribution. 
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The effect of the kinematical smearing on p 2 may be seen by comparing the distri- 
butions for the generated events (middle row) to the ones for measured events (bottom 
raw) where the acceptance and smearing were included. The smearing of pi increases with 
increasing pi; it is about 0.006 GeV 2 for the coherent samples (pf < 0.02 GeV 2 ) and about 
0.03 GeV 2 for the incoherent samples (pi > 0.1 GeV 2 ). 

The analysis of the p decay distributions allows us to study spin-dependent properties 
of the production process [32], in particular the polarization of p. Usually the p° decay 
angular distribution W(cos8, <p) is studied in the s-channel helicity frame, which is the 
most convenient for describing the p decay after photo- and electroproduction [36] . The p° 
direction in the virtual photon-nucleon centre-of-mass system is taken as the quantization 
axis. The angle 6 is the polar angle and <f the azimuthal angle of the tt + in the p° centre- 
of-mass system. The cos 9 distributions for pions from p° decays are shown in Fig. 4 for 
the lead target. The distributions for longitudinally (dashed lines) and transversely (dotted 
lines) polarized parent p 0, s are markedly different. Their sum is also indicated. 

Usually fits to the combined cos 6 distributions are performed in order to determine 
the density matrix element t[}q (cf. Eq. 11), which can be identified as the probability that 
the p° was longitudinally polarized. For exclusive p° production the approximate s-channel 
helicity conservation (SCHC) is observed [37, 18], i.e. the helicity of p° is predominantly 
equal to that of the virtual photon. Assuming SCHC and using the fitted Tqq one can 
estimate the ratio R = o~l/°~t for exclusive virtual photoproduction ([32]). In COMPASS 
we plan to measure R as well, and study its Q 2 - and ^4-dependence, which is expected to 
reflect the strength of nuclear absorption. 

As Q 2 increases, the approach to the CT limit is expected to be different for VMP 
by longitudinally polarized virtual photons from that by transversely polarized photons. 
Therefore, in order to increase the sensitivity of the search for CT, we propose another 
method. It consists in studying ^-dependence of the cross sections for samples with different 
p° polarizations, which will be selected by cuts on cos 9. For instance, after applying the 
cut | cos# |> 0.7 the fraction of accepted events with p\ is (80-95)% depending on the 
simulation, whereas for the cut | cos 8 |< 0.4 the fraction of events with pj, is (75-92)%. 
For an approximate SCHC, such cuts will allow us to select the samples with enhanced 
contributions of the events initiated by the virtual photons of a desired polarization. Studies 
of the samples with different polarizations of the virtual photons are important for the clear 
demonstration of CT. 

Another aspect which is important for CT studies, is the covered range of the co- 
herence length l c (cf. Section 1). In Fig. 5 we present the plot of l c vs. Q 2 for a sample 
of accepted events. The effects of initial and final state interactions in the nuclei vary at 
small l c values [15]. Therefore, it was suggested that in order to disentangle effects due to 
CT from those caused by the modified absorption at small l c values, one should study A- 
and Q 2 -dependences of cross sections at fixed values of l c . This approach will be possible, 
if large statistics data were available. 

For a limited statistics, a possible solution to avoid the mentioned effects is to use 
the combined data at l c values exceeding the sizes if the target nuclei. The radius of the 
carbon nucleus is < > 1 / 2 = 2.5 fm and that of the lead nucleus is < 7p b > l / 2 = 5.5 fm 
[38]. Therefore, one may use the selection l c > lf in ~ 2- < r£ b > l / 2 = 11 fm. The value of 
l™ m is indicated in Fig. 5 by arrows. About a half of events survive the cut on l c . These 
events cover the range of Q 2 < 6 GeV 2 , which is expected to be sufficient to observe CT. 

The estimated values and statistical precision of R T , the ratio of the nuclear trans- 
parencies for lead and carbon, are presented for different Q 2 bins in Fig. 6 for pi < 0.02GeV 2 
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Table 3: Q 2 bins used for the determination of R T . 



Kill VIllTYlhDT 1 

J_>111 IlUIilUt:! 


D 2 Kin 
\cg Dill 


<r n 2 -> 






[GeV 2 ] 


[GeV 2 ] 




1 


2-3 


2.4 


0.016 


2 


3-4 


3.4 


0.022 


3 


4-6 


4.8 


0.031 


4 


6-9 


7.2 


0.048 


5 


9-12 


10.2 


0.072 


6 


12-20 


14.8 


0.11 



and in Fig. 7 for p 2 > 0.1 GeV 2 . The Q 2 bins are specified in Table 3. Each figure comprises 
predictions for two models, CT and NA, and for 6 different samples of accepted events for 
each model. For each sample a set of "measurements" in different Q 2 bins is shown. Sets 
A and B were obtained using the standard selections for the MT and FT triggers, respec- 
tively. Four remaining sets were obtained for the FT trigger with additional selections: C 
with | cos | < 0.4, D with | cos 6* | > 0.7, E with | cos 6* | < 0.4 and l c > 11 fm, and F with 
| cos 0\> 0.7 and l c > 11 fm. Note that for sets E and F only three lower Q 2 bins appear (cf. 
Fig. 5). One expects large differences in i? T for the two considered models. For coherent 
samples i? T ~ 5 for CT model and ~ 1 for NA model. At Q 2 ~ 5 GeV 2 the precision of 
the measurement of R T for coherent events will be better than 17%, even for the restricted 
samples E and F, thus allowing excellent discrimination between the two models for the 
nuclear absorption. For the incoherent events the power to discriminate models by R T 
measurements will be more limited. 

5 Comparison with previous experiments 

In this section we compare COMPASS capabilities to demonstrate CT to those of 
previous experiments in which exclusive p° leptoproduction on nuclear targets was studied. 
We concentrate on this reaction, as among different exclusive VMD channels it has the 
largest cross section. For the comparison we have selected experiments which covered Q 2 
range extending to large values, bigger than 2 GeV 2 . This condition is satisfied for the 
following experiments: HERMES [15], NMC [25] and E665 [24]. 

The published results from the HERMES experiment concern the incoherent exclusive 
p° production on 1 H, 2 H, 3 He and 14 N targets. The electron beam energy was 27.5 GeV 
and the covered Q 2 and l c ranges are 0.4 < Q 2 < 5.5 GeV 2 and 0.6 < l c < 8 fm. Like in all 
fixed-target experiments with lepton beams, the covered kinematic range is such that for 
a given Q 2 the average value of l c decreases with increasing Q 2 . The highest values of the 
HERMES data where one may expect to observe the onset of CT are correlated with the 
small values of l c below 2 fm. The increase of the nuclear transparency observed for 14 N is 
explained as due to the reduced coherence length, without resorting to CT. 



14 



Wider kinematic ranges were coverd by experiments using high-energy muon beams. 
NMC has published the data on coherent and incoherent exclusive p° production on 2 H, 
C and Ca targets. The muon beam energy was 200 GeV and the data cover the ranges 
2 < Q 2 < 25 GeV 2 and 1 < l c < 30 fm. For the NMC kinematic range, the large values of 
l c which exceed the diameter of the largest target nucleus (Ca), l c > 7 fm, correspond to 
Q 2 < 5.5 GeV 2 . In principle, in this range it is possible to study Q 2 dependence of nuclear 
absorption in order to observe CT not obscured by effects of small l c . However, due to the 
moderate statistics of the data, such detailed analysis was not possible for NMC. 

The most favourable kinematical conditions were realised in the E665 experiment due 
to the high muon beam energy of 470 GeV. The experiment has published the data on 
incoherent exclusive p° production on 1 H, 2 H, C, Ca and Pb targets. The data were taken 
in a wide Q 2 range, including also very small values Q 2 > 0.1 GeV 2 , and in a wide range 
of l c , 1 < l c < 200 fm. At small and moderate Q 2 values, < 3 GeV 2 , which correlate with 
l c values exceeding the diameter of the lead nucleus, precise measurements of the nuclear 
transparency were obtained. They indicate strong nuclear absorption. At the highest Q 2 
bin, Q 2 > 3 GeV 2 , the nuclear transparency increases in qualitative agreement with CT. 
However, the data in this bin correspond to a wide range of l c , 1 < l c < 40 fm, and therefore 
could be affected by effects due to the reduced coherence length. The low statistics in the 
large Q 2 bin did not allow more detailed studies. 

The COMPASS data at medium and large Q 2 (FT trigger) will cover the kinematic 
range similar to that of the NMC data. The expected statistics for the carbon target will 
be about 2 orders of magnitude higher than that of the NMC data for the same target. 
This increase is in particular due to higher beam intensity and larger acceptance in the 
COMPASS experiment. Similarly like in the E665 experiment COMPASS will use also 
lead target and will extend measurements to small Q 2 . The statistics of COMPASS data 
will be significantly higher also than that of E665. For example, for incoherent events at 
Q 2 > 3 GeV 2 the ratio of nuclear transparencies for lead and carbon will be determined with 
accuracy better than 3%, comparing to about 30% for E665. Another more sensitive method 
to study nuclear effects in exclusive p° production, the one by using coherent production, 
will be employed in COMPASS as well. Due to large statistics splitting of COMPASS data 
in several Q 2 and l c bins as well as the selection of events with longitudinal or transverese 
p° polarization will be possible. 

6 Open questions 

In addition to the p° production, we will study exclusive production of J/ip, 4>, ip' and p' 
and their ratios in order to demonstrate CT. For these channels the analysis will be similar 
to that for p°, but the rates will be lower due to the smaller production cross sections. The 
estimates of the required period of data taking and of the target thicknesses are subject of 
further analysis. 

We consider also a possibility to study CT via coherent production of hadron pairs 
with large relative transverse momenta. The experimental procedure would be similar to 
that used in the Fermilab experiment E791 [26], in which CT was observed in coherent 
diffraction of pions into high-p^ di-jets. We will select a sample of events with high-pt 
hadron pair and an additional requirement that all observed charged hadrons carry at least 
90% of the energy of 7*. An estimate of t will be obtained using momenta of all measured 
charged hadron tracks. The expected resolution in t, which will be crucial for selection of 
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the coherent sample, will be determined from a dedicated simulation and reconstruction of 
events with high-p t hadron pairs. 

Finally we mention few experimental aspects which have to be further investigated. 

Here we propose to position the nuclear targets downstream of the polarized target. 
Then the CT physics could complement the standard COMPASS program by using the 
downstream nuclear targets simultaneously with the polarized target. Similar target con- 
figuration was already used during the 1985 running of the EMC experiment, when the 
polarized target was used together with several nuclear targets situated downstream of PT. 

The effect of the downstream nuclear targets on the hadrons coming from the polarized 
target should be small. For instance, for the reconstructible tracks of hadrons from PT, only 
about 4% of slow pions from D* decays and about 1% of hadrons from D° decays will be in 
the acceptance of the proposed nuclear target [39]. For the scattered muon coming from PT 
and passing through the nuclear target, the additional multiple Coulomb scattering should 
result in a small increase of the errors of the reconstructed position of the event vertex. 

Assuming a dedicated run with unpolarized targets, different target configurations 
for CT studies are possible. With more available space one could envisage targets at least 
about 4 times thicker than the thin targets presented earlier. Thus the four-fold reduction 
of the data taking time will be possible. With an even larger amount of available beam 
time, one could expose more nuclear targets of different materials, and in addition the liquid 
hydrogen or deuterium targets. The latter ones are also important for measurements of off- 
forward parton distributions (also referred to as skewed parton distributions or generalized 
parton distributions) [40]. For a dedicated run with unpolarized targets, one could take 
advantage of a possibility of adding a recoil detector. This is crucial for certain processes, 
e.g. Deep Virtual Compton Scattering (DVCS), but in any case it will help to reduce the 
non-exclusive background for any exclusive process. 

7 Conclusions 

We have simulated and analysed exclusive p° muoproduction (jxA — > n p°A and \xA — > 
p p°N(A — 1)) at the COMPASS experiment using thin nuclear targets of carbon and lead. 
For muon beam energy of 190 GeV and a trigger for medium and large Q 2 , the covered 
kinematic range is 2 < Q 2 < 20 GeV 2 and 35 < v < 170 GeV. 

Good resolutions in Q 2 , l cl t (p 2 ) and cos 8 are feasible. An efficient selection of 
coherent or incoherent events is possible by applying cuts on p 2 . In order to obtain the 
samples of events initiated with a probability of about 80% by either 7£ or 7^, the cuts 
on the p° decay angular distribution of cos 8 will be used. The search for CT could be 
facilitated by using the events with l c values exceeding the sizes of the target nuclei. The 
fraction of such events is substantial and the covered Q 2 range seems sufficient to observe 
CT. 

We showed high sensitivity of the measured ratio i? T of nuclear transparencies for 
lead and carbon for different models of nuclear absorption. Good statistical accuracy of the 
measured R T may be achieved already for one year of data taking. These measurements, 
taken at different Q 2 intervals, may allow to discriminate between different mechanisms of 
the interaction of the hadronic components of the virtual photon with the nucleus. 

In conclusion, the planned comprehensive studies of exclusive vector meson production 
on different nuclear targets at the COMPASS experiment should unambiguously demon- 
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strate CT. 
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Figure 1: Diagram of exclusive meson muoproduction \i N — ► fi N V , where meson V is the 
only produced particle. Corresponding four-momenta are indicated in brackets (cf. Table 
1) • 
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Figure 2: Inelasticity distribution. The SMC preliminary results [29] for jiN — > fip°N . 
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Figure 3: Distributions of p\ for the lead (left) and carbon (right) targets. The distributions 
a) and b) at the top are for the generated events and the wide range of pi , whereas those 
in the middle, c) and d), correspond to the range of low p\. The bottom distributions, e) 
and f), are for the accepted events, with the kinematical smearing taken into account. The 
arrows show the cuts described in the text. 
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Figure 4: Distributions of cos 6 for the pions (tt + ) from the decays of p° mesons produced 
on lead, for generated (left) and accepted (right) events. For the latter ones the effects of 
the kinematical smearing were taken into account. The dashed- and dotted-line histograms 
are for the longitudinally and transversely polarized p°'s, respectively, and the solid-line 
histograms are for their sum. 



23 




Figure 5: Correlation of l c and Q 2 . The plot a) and its projections b) and c) were obtained 
for the accepted events from p° production on lead. The arrows show the cut discussed in 
the text. 
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Figure 6: Expected statistical precision of the measurement of the ratio i? T of nuclear 
transparencies for lead and carbon for pi < 0.02 GeV 2 and for different Q 2 bins, for CT 
(upper band) and for NA (lower band) models. The Q 2 bins are defined in Table 3. Sets 
A, B, C, D, E and F correspond to the conditions specified in the table below the plot. 
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Figure 7: Expected statistical precision of the measurement of the ratio R T of nuclear 
transparencies for lead and carbon for p 2 > 0.1 GeV 2 and for different Q 2 bins, for CT 
(upper band) and for NA (lower band) models. The Q 2 bins are defined in Table 3. Sets 
A, B, C, D, E and F correspond to the conditions specified in the table below the plot. 
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